Detect Spinons via Spin Transport 
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Existence of spinons is the defining property of quantum spin Uquids. These exotic excitations 
have (fractionalized) spin quantum number and no electric charge, and have been proposed to form 
Fermi surfaces in the recently discovered organic spin liquid candidates. However direct probes for 
them are still lacking. In this paper we propose to experimentally identify the spinons by measuring 
the spin current flowing through the spin liquid candidate materials, which would be a direct test for 
the existence of spin-carrying mobile excitations. By the nonequilibrium Green function technique 
we evaluate the spin current through the interface between a Mott insulator and a metal under a 
spin bias, and flnd that different kinds of Mott insulators, including quantum spin liquids, can be 
distinguished by different relations between the spin bias and spin current. In the end we will also 
discuss relations to experiments and estimate experimentally relevant parameters. 

PACS numbers: 75.10. Jm,72.25.Mk,73.23.-b, 73.40.Rw 



Quantum spin liquid(QSL) was first proposed by An- 
derson as an alternative ground state against long range 
magnetic order in frustrated magnets^. In these systems 
competing spin exchange interactions result in a large 
degeneracy of classical ground states, and quantum fluc- 
tuations among these states destroy long range symme- 
try breaking order A particular kind of quantum spin 
liquid, the resonant valence bond(RVB) state, has also 
been proposed to be the key to the high-temperature su- 
perconductivity in the cuprate materials <^After decades 
of intense research many numerical evidences of QSL 
ground states have been found in semi-realistic lattice 
models^"—, and many artifical parent Hamiltonians for 
spin liquids have been constructedi^"— . 

On the other hand, the experimental realization of 
spin liquids in more than one spatial dimensions re- 
mains challenging until several candidate materials have 
been discovered recentlyi^ii^ Two two-dimensional(2D) 
triangle lattice organic salts EtMe3Sb[Pd(dmit)2]2 and 
K-(BEDT-TTF)2Cu2(CN)3, the kagome lattice herbert- 
smithite ZnCu3(OH)6Cl2 and a three-dimensional hyper- 
kagome lattice Na4lr308 are found to be the most 
promising candidates of QSL<^"— Despite of structural 
distinction, they are all Mott insulators with compet- 
ing interactions and show no magnetic order down to 
much lower temperature than their exchange interaction. 
Current measurements of magnetic susceptibility, spe- 
cific heat, thermal transport and neutron scattering have 
provided vital information about the properties of these 
materials j^^— However, there is still no definitive exper- 
iment for the identification of quantum spin liquidsi^ 

One of the most significant features of quantum spin 
liquids is that they have exotic excitations called spinons 
which are uncharged and usually spin- 1/2 mobile parti- 
cles, may obey bosonic or fermionic statistics and may 
have a gap or not^^ The fermionic spinons may form 
Fermi surfaces and are generally accompanied by an 
emergent gauge field i^^'^^ This "spinon Fermi sea" state 
has received strong support from the observations of 



metallic-like specific heat and thermal conductivity^ in 
the organic candidates at low temperatures. However 
these experiments do not provide a direct proof that the 
mobile and possibly fermionic low energy excitations are 
spinons. A more reliable proof for the spinon Fermi sea 
would be the metallic-like spin transport in these Mott 
insulators. 
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FIG. 1: (Color online). Schematic plot of a four-terminal 
device for the detection of spin current through a spin liquid. 
Four small blue bars indicate ferromagnetic electrodes. A 
spin polarized current is injected into the right lead(R-lead) 
by a current source(/e), thus creates a spin bias and drives 
a spin current through the spin liquid middle region. The 
spin bias induced by the spin current in the left lead(L-lead) 
is then measured by a voltmeter(1/). We will also consider 
other kinds of Mott insulators as the middel region instead of 
spin liquids. 



In this paper, we propose a four-terminal measurement 
of spin current through a spin liquid material (Mott in- 
sulator) as the evidence for the existence of spinons. The 
proposed four-terminal device consists of a spin liquid 
material coupled to the left and right normal metal leads 
and each lead couples to two ferromagnetic(FM) elec- 
trodes as shown in Fig. [1}^ A current source is added 
between the two right FM electrodes. This is used to 
create a spin polarized current flowing from one right 
FM electrode through the right lead to another right 
FM electrode, leading to a spin-solved chemical poten- 
tial (i.e. a spin bias Vr) in the right lead. Then this 
spin bias will drive a spin current flowing from the right 
lead through the spin liquid and finally into the left lead 
by the spinon-electron spin-exchange interaction at the 
interfaces between the spin liquid and the leads. At last. 
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a voltmeter is connected to the two left FM electrodes 
which are in contact to the left lead. The voltmeter is 
used to measure the spin bias Vl created by the spin cur- 
rent in the left lead. Here, the spin bias Va (a = L, R) 
is defined as the difference between the spin-f chemical 
potential and the spin-^ chemical potential fiai in 
the a-lead, i.e. Va = Mat ~ fJ-ai-^^^ 

In the rest of this paper, we will first establish the gen- 
eral result of the spin current through the interface be- 
tween a Mott insulator in the middle region and a metal- 
lic right lead with a spin bias (see Fig. [1]) by the nonequi- 
librium Green function technique. We will then apply the 
general formalism to show that different Mott insulators 
can be distinguished by different relations between the 
spin current and the spin bias as well as temperature. In 
the end we will also discuss relations to experiments and 
estimate experimentally relevant parameters. 

The model Hamiltonian and formulation. In our theo- 
retic analysis, we consider the model of a Mott insulator 
(a spin liquid or a coUinear antiferromagnet) coupled to 
two normal leads under a spin bias on the right lead. The 
Hamiltonian of the system is given hj H = Hq+Hm+Hj, 
where Hq, Hm, and Hj are the Hamiltonians of the leads 
(Metal), the middle region (Mott insulator), and the in- 
terfaces respectively. Hq and Hj are assumed to be 



the right lead is. 



(1) 



and 



(2) 



where 



with SM{f) the (dimensionless) spin operator of the mid- 
dle region at the position r. c^^^{c\^^) is creation (an- 
nihilation) of the spin-cr electron in the Q;=(L,R)-lcad. 

— fJ'aa, where Sak is the electron dispersion 
relation and A''acr=(t,i) the spin dependent chemical 
potential in the a-lead. The spin exchange interaction 
constants J/ is determined by the interface properties of 
Mott insulator and metal<^ ™^6ans integral over the 
interface. Hm depends on the type of the Mott insulator 
we consider and will be specified later. We emphasize 
that there is no single electron tunneling term in Hj be- 
cause the middle region is a Mott insulator. 

Due to the spin exchange interaction in Hj, the spin 
current can flow from the normal lead to Mott insulator 
and vice versa. When the right lead is under a spin bias, 
the spin current Is flowing into the middle region from 



with 



Is = - <-S^>=t<[S^,H{t)]> 
at 

= iJiY.< [S^,SR{ro) ■ SMiro)] > 

ro 

= J/^i?e(r<(ro,t,t)), 

''0 



(3) 



{h/2){N^-N^) = [hmY.ki^'Rk^c^.k^ - 

Here, we have used the fact that 
and defined r<(ro,t,t) = i < 
SM{ro,t)c\^{rQ,t^}cj^^{r(i,t) > with S]^j{ro) = 

SUro)-iSl,{ro). 

In order to solve the Keldysh Green func- 
tions above, we first apply equation of motion 
technique to solve r*(ro, fe, fc', r, r') = —i < 

r,{5,-,(ro,r)c^,,(r)4,,,(r')} 

By keeping the lowest order terms of J/ we have 

r*(ro,fc,fc',T,T') 



-iJi 



2N„ 



xg*j^-f{k' ,ti,t') exp[-i{k - k') ■ ro], 

where x\ro,r'Q,T,Ti) ^ -i < Tc[SJi{ro,T)Sji{r'„,Ti)] > 

and 5*j,(fc',T,T') = < rc[c^;,,,(T)4^,^(T')]^ > are 
contour-order Green functions for spin operator in mid- 
dle region and free electrons in right lead, respectively. 
T<{ro,k,k',t,t') = i < 5^/ (^"0, Ocflfe'T(*')cj?fci(*) > can 
be obtained by analytic continuation from F*, and a 
Fourier transform then produces F^(ro,t,i). Plug the 
result into Eq. (jS]) we have 



^■■^ = AN^ /V2 E ^M{qARk't-^B.kl + V)Sq^+k^ 



q.k.k' 



-k\ 



X { [1 + ns (Cijfc't - iRki + V)]nF {-£,Rki)nF i^Rk't) 
-nB {£,Rk't~ ^Rki + V)nF {^Rki)nF (-Ci?fc't ) } (4) 

where Nm and Nr are the number of unit cells in the 
middle region and right lead, respectively. Here N± 
is the number of transverse mode (parallel to the in- 
terface) and V = fiR-f ~ /iijj, is the spin bias in the 
right lead. nsiLo) and nF((,Rka) arc the Bose and 
Fermi distribution functions, respectively, bq^^k^-k' 
indicates transverse(parallel to the interface) momen- 
tum conservation. The power spectrums Am are de- 
fined as ^A^(q,a;) = J dt < S^j{~q,t)Sjj{q,0) > 
exp{iu!t)/nBi(jj), where S^j{q,t) is the Fourier transform 
ofS^,{r,t) = SfAr,t)±Sl{r,t). 

We note that the behavior of the spin current is mainly 
determined by the power spectrum Am of the mid- 
dle region. Under appropriate conditions the momen- 
tum integrations over q, k, k' can be approximately sep- 
arated, and the transverse momentum conservation fac- 
tor will provide only a constant factor— . The 
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Fermi (Bose) function will be treated exactly at zero tem- 
perature and expanded in series of Vj (ksT) at finite tem- 
perature. In the following parts of this paper we will ap- 
ply Eq. (HI and analyze several kinds of Mott insulators 
as the middle region, including several spin liquids. 

Spin liquids. Spinous in spin liquids may or may not be 
gapped. The gapped spin liquids will have exponentially 
vanishing spin transport at low temperature and small 
spin bias. We therefore restrict ourselves to the type of 
QSLs with gapless fermionic spinous. We describe such 
spin liquids by the following mean-field Hamiltonian 

k,(T 

where / are fermionic spinous and Cfc = Cfc — fJ-s with 
fis the spinon chemical potential. The spinon dispersion 
ek may have a Fermi surface (the "spinon Fermi sea" 
stat o^^'^^ ) or Dirac points at Fermi level^. 

For illustration purpose we first consider a one- 
dimensional(lD) spinon Fermi sea state. In ID case 
the transverse momentum conservation factor in Eq. Q 
does not exist, and the momentum integrations can 
be done separately. At T = and >> w, since 
Cq = h^q^ / {2ms) with the spinon effective mass, the 
power spectrum Am divided by N]^ and integrated over 
momentum q (the "density of states of spin excitations" ) 
is 

R q 

Replace it in Eq|31 we find spin current oc V^'^ at T = 0, 
and Is oc {ksTfY at T > with V « ksT. 

Now we come to the two dimensional spinon Fermi sea 
case, which is the most relevant to the 2D organic spin 
liquid candidate materials. 

The spinon dispersion is = h'^q^ / {2ms) with 
the spinon effective mass, the density of states of spin 
excitations is 

^^AA/(g,w) = 2nN''s{E'p)u^oj. 

R q 

Thus the spin current /g oc F'^ at T = (see Fig 12]), and 
Is oc {kBT)'^V at T > with V « kgT (see FigE]). 

In the two dimensional Dirac spin liquid case^^, the 
low energy spinon dispersion is eg = ±hvp\q — qp\^ with 
vp the Fermi velocity and qp the Fermi vector. 

At T=0, since Ylq ^m('?, w)/iV|. oc lo'^ , the spin current 
Is oc V^. At T > 0, V « ksT, expanding Is in the 
scries of V/{kBT) directly, we find /, oc {kBT)'^V. 

Collinear antiferromagnetic Neel order. Antiferromag- 
netic(AFM) order is a common competitor for quantum 
spin liquids. Let us now consider the simplest AFM or- 
dered state, the collinear AFM Neel order on a bipartite 
lattice, and show that it has different spin transport be- 
havior compared to the previous spin liquid cases. We 
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FIG. 2: (Color online) In 2D spinon Fermi sea case, the scat- 
tered dots are the spin bias V dependence of Is/V^ at zero 
temperature (fitted by simple line). The cut point at origin 
directly indicates 7s oc 1/'^. 
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FIG. 3: (Color online) In 2D spinon Fermi sea case, (a) the 
scattered dots are temperature T dependence of la/T {V the 
spin bias on the right lead). They can be fitted by simple 
lines(y oc x). (b) The spin current Is versus spin bias V on 
the right lead at difi'erent temperatures T (fitted by simple 
lines). The cut point at origin indicates Is oc V. (a) and (b) 
together indicate Is oc T^V at high temperature ksT >> V. 

describe the spin excitations by a linearized spin wave 
Hamiltonian, 

Hm = Eo{ala^ + blb^ + Jka^b_^ + -f^albl^) 

k 

where Eq ~ 2ZS\J\ and 7^. ~ X]<sCos(fc • 6)/Z with 5 
sums over the nearest-neighbor lattice sites. Here Z is the 
coordination number. S is spin quantum number and J is 
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TABLE I: Behavior of the spin current through the interface 
between different Mott insulators(rows) and a metallic lead 
with respect to the spin bias V and temperature T. 





T = 


kgT » V 


1-d spinon Fermi sea 


oc V -" 




2-d spinon Fermi sea 


(X V 


oc {ksTYV 


2-d Dirac spin liquid 


cx V 


oc {kBTYV 


coUinear antiferromagnet 


oc V\V\^ 


oc {kBTyv 



TABLE IL Numerical estimates of the induced spin bias in 
the left lead with different kinds of Mott insulator middle 
regions(rows) at three different temperatures T = OK, IK, 
and IQK. Other parameters used are given in the main text. 
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the spin exchange constant. 6^ {b\.) and (1^), are the 
annihilation (creation) of Holstein-Primakoff bosons^ on 
B-sublattice and A-sublattice, respectively. 

Since spin wave dispersions is Eq = 4'\/2>5'| J|ga with 
a the lattice constant on square lattice, the density of 
states of spin excitations is 

R q 

The spin current 1^ oc V\V\^ and h oc {ksTfY at T=0, 
and T > with V « ksT, respectively. 

All the cases we have discussed above are summarized 
in Table n 



Numerical estimates of the spin current. We use the 
following estimates of the parameters^, with the inter- 
face exchange coupling Jj = lOmeV and spin bias V = 
IkB ~ O.lmeV. The metallic conductivity p = 10~*J7-m, 
lattice constant a = 0.3nm, and Fermi level Ep = leV. 
The effective spinon mass w lOm^ and Fermi level 
Ep = lOmeV. The bias induced in the left lead is evalu- 
ated by Vl ~ {2e/h){Is/N±) * p*ls with N± the number 
of transverse mode (see Table . 

Discussions. Many factors ignored by our analysis may 
affect the results. First we have assumed the conservation 



of the z-component of spin in the entire system, so the 
spin current is well-defined. However in the real materi- 
als spin-orbit coupling(SOC) will generically be presenli^. 
We hope our results can still be applied to such systems 
if the linear dimensions of the sample are much smaller 
than the inverse of SOC. For the same reason we did not 
considered non-collinear AFM orders, e.g. the 120° order 
on triangular lattice. 

Secondly we have ignored the emergent U{1) gauge 
field in the spinon Fermi sea and Dirac spin liquid cases. 
It is well-known^i^ that coupling to this gauge field 
can significantly change the low energy behaviors of 
the (spinon) Fermi sea. However such effects have not 
been found in the specific heat and thermal conductiv- 
ity measurements of the organic spin liquid candidates, 
the conventional Fermi liquid behaviors were observed 
instead2ii^i^. We therefore believe our results are still 
valid in these materials. The effect of U{1) gauge field 
is an interesting theoretical question and will be left for 
future studies. 

Finally we have assumed a clean and free spinon or 
magnon system in the middle region, without any scatter- 
ing of spinous or magnons by interactions among them- 
selves or impurites. We think this is not a serious problem 
for experiments, according to the large value of Ipm of 
the experimentally estimated spinon mean free patb^S. 

In summary, we propose to experimentally identify the 
spinous by measuring the spin current flowing through 
the spin liquid candidate materials, which would be a di- 
rect test for the existence of spin-carrying mobile excita- 
tions. By the noncquilibrium Green function technique 
we evaluate the spin current through the interface be- 
tween a Mott insulator and a metal under a spin bias. It 
is found that different kinds of Mott insulators, including 
quantum spin liquids, can be distinguished by different 
relations between the spin current spin bias as well as 
temperature. We hope our results can stimulate more 
experimental studies of the spin liquid candidate mate- 
rials and further promote the exchange of ideas between 
different fields (e.g. spintronics and strongly correlated 
electrons) in condensed matter physics. 

This work was financially supported by MOST 
of China (2012CB921303, 2009CB929100 and 
2012CB821402), NSF-China under Grants Nos. 
11074174, 11121063, 91221302 and 11274364. 
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I. SPIN CURRENT 

f>r^ ' In the main text, we have introduced 
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-nsi^Rk't - ^Rki + V)nF{^Rki)nF{-^Rk't)} (1) 



where Nm and Nn is the number of unit ceUs in the middle region and right lead, respectively. Here N± is the 
number of transverse mode(parallel to the interface) and V = fiR^ — fiRi is the spin bias in the right lead. nB{uj) and 
^ "-fCC-R/cct) are the Bose and Fermi distribution function, respectively. (5qj^+fcj^-fc'^ indicates transverse (parallel to the 
^ • interface) momentum conservation. The power spectrums Am are defined as Am{(1, i^) = J dt < S'^j(— q, t)S~^.i{q, 0) > 
^ ' exp(ia;t)/nB(aj), where S^^{q,t) is the Fourier form of S^j{r,t) = S'^.j{r,t) ± S^j{r,t). 

In the main text, we have assumed that the transverse momentum conservation term can be separated properly. 
This will be explained it in detail in the rest of this material. 



^ , II. SPIN LIQUID 

o , 

^ We describe gapless spin liquids by the following mean-field Hamiltonian 



k,(T 

where / are fermionic spinons and Ck = ^k — fJ-s with fis the spinon chemical potential. In this case, Sa/ (r) = 
i ^ o'fj^fj^i fli{''')f^,{r), with cr the Pauli matrix and /i =t,i- The spinon dispersion may have a Fermi surface (the 

"spinon Fermi sea" stateiiS) or Dirac points at Fermi level^. 
Because 



we get 



y^_^AM{q,uj)Sq^+k^-k', 

E nF{C,p-q)nF{-Cp)6{uj - Cp + Cp-q) 



E ^P±-p'^+k±-k\ 



nF{Cp')nF{-Cp)S{^ ~ Cp + Cp') 



p ,p 
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Now we come to the two dimensional spinon Fermi sea case, which is the most relevant to the 2D organic spin liquid 
candidate materials. If the energy of spinon excitations (^k and spin-f {],) electron excitations ^RkfH) a-re much smaller 
spinon Fermi level fis and spin-f (i) electron Fermi level /x_R,-f(4.) i 

p± — + k± — k'j_ = p cos 01 — p' cos d'l + k cos 62 — k' cos d'2 

~ kp cos 61 — kp cos O'l + cos O2 — kp^ cos d'2 . 

In the two dimensional Dirac spin liquid case^, the low energy spinon dispersion is = ±hvp\q — qfI, with vp the 
Fermi velocity and qp the Fermi vector, 

p± — p'^ + k± — k'j_ = p cos 61 — p' cos 6'i + k cos 02 — k' cos 6*2 
« k'j^^ cos 62 — kp^ cos d'2 ■ 
Since p± — p'^ + k± — k'j_ is the only angular dependent term, we can express the spin current as 

" 4jV2jv2 ^A/(g,4flfe'T - ^Rki + V) 

x{[l + nBi^Bk'f - £,Rki + V)]np{-^iiki)nF{^Rk't) 

-nB{£,R.k't - ^RH + V)nF{^Rki)nF{-£.Rk't)}, (2) 
with Qi — J ddid9'^d92dd'2{cosdi — cos9'^ + kp^/kp cos02 — kp^/kp cos 6'2) / {2tt)^ for spinon Fermi sea case and Qi = 
J dd2dd'2{cos92 — kp^/ kpi^ cos d'2) /{2tt)'^ for U(l) Dirac spin liquid, respectively. 
Thus the transverse momentum conservation term is separated independently. 
We note that in one dimension case, Qi = I and the number transverse mode N± = 1. 

III. COLLINEAR ANTIFERROMAGNETIC NEEL ORDER. 

Let us now consider the simplest AFM ordered state, the coUinear AFM Neel order on a bipartite lattice, and 
show that it has different spin transport behavior compared to the previous spin liquid cases. We describe the spin 
excitations by a linearized spin wave Hamiltonian, 

k 

where Eq = 2ZS\J\ and 7fc = -^^^co^ik ■ 5) with 5 sums over the nearest-neighbor lattice sites. Here Z is the 
coordination number, S is spin quantum number and J is the spin exchange constant. In this case, S~^{r) = 
SUr)+iSl{r)^a^r) + b{r) , 5,-,(r) = Sl,{r) - tS^ir) = a{r) + bHr) and Slj{r) = a\r) a (r) - bHr)b{r), where 
a(r) [a'l'(r)] and 6(r) [fo^('r)] are the Fourier form of a{k) [a^(A;)] and b{k) [b''{k)], respectively. Here, b/, (6^^) and 
(a^) are the annihilation (creation) of Holstein-Primakoff bosons'* on B-sublattice and A-sublattice, respectively. 
If the wave vector of spin wave is much smaller than spin-f (J,) electron Fermi vector, i.e. ki << kp^{kp^), 

q± — k± + k'j_ = q cos 9i — k cos 02 + k' cos ^3 « kp^ cos 92 — kp^ cos 6*3. 

The spin current is 

R q,k,k' 

x{[l + nsi^Rk't - fflfci + V)]np{-(,Rki)np{(,Rk'f) 
-nsi^Rk't - ^Rki + V)nF{^Rki)nF{-£,Rk't)}, (3) 
with = / d62d9r^6 {cos 92 — k p^ / kp^ cos 9^) /{2Tr)'^. Thus the transverse momentum conservation term is separated. 
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